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Abstract Isospin-dependent Quantum Molecular Dynamics model (IQMD) has been applied to investigate the 


Pygmy Dipole Resonance (PDR) and Giant Dipole Resonance (GDR) in Ni isotopes by Coulomb excitation. By 


Gaussian fitting to the photon emission spectra, the peak energies and strengths of PDR and GDR are extracted. Their 


sensitivities to impact parameter, incident energy and the symmetry energy are discussed. By the comparison of 


energy-weighted sum rule (EWSR) with the data and other calculations for “Ni, the parameters of density-dependence 


of symmetry energy in the IQMD are constrained. In addition, the N/Z dependence of PDR and GDR parameters of Ni 


isotopes are investigated, and the results that the EWSR increases linearly with the N/Z are obtained. 
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1 Introduction 


In the past decades, many works have been done on 
Pygmy Dipole Resonance (PDR) and Giant Dipole 
(GDR) both 
experimentally“). PDR or GDR is classic collective 


Resonance theoretically and 
mode of nuclei, which can be considered as the 
vibration of valence neutrons against the nuclear core 
or neutrons against protons, respectively. 

The results of many works show that the PDR 
plays an important role in neutron-capture rates in the 
r-process, nucleosynthesis, the radiative neutron- 
capture cross section on neutron-rich nuclei, as well as 
the photodisintegration of ultra-high energy cosmic 
rays”, Recently, the correlation between PDR and 
symmetry energy has also been studied!!?"9). 

In the present work, we try to apply the 
Isospin-dependent Quantum Molecular Dynamics 
(IQMD) model to investigate PDR and GDR in Ni 


isotopes. 


2 Model and methods 


In this section, we introduce the model and formulae 
which we used in the present calculation as well as the 
method how we choose the valence neutrons. 


2.1 Isospin-dependent quantum molecular 


dynamical model 


The Isospin-dependent Quantum Molecular Dynamics 
model which is based on QMD model is a kind of 
Monte-Carlo erent: Tt 
extensively applied in heavy-ion collision dynamics. 

In our IQMD model, the mean field can be 
written as follows: 


transport mode has been 


U(p)=U™ you ye us y™ (1) 


Here, USY, US, UX, U%™ and UM is the 


Skyrme potential, Coulomb potential, Yukawa 


potential, symmetry potential and Momentum 


Dependent Interaction (MDI), respectively. 
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The Skyrme potential is given as follows: 
U =al p/P) + B(P/ Po)” (2) 


Where pois the saturation nuclear density (po=0.16/fm*) 
and p is the nuclear density. Parameters [a, p, o] which 


Table 1 Parameters [a, /, o] for different EOS 


represent the equation of state are given in Table 1. In 
the table, S(M) represents the soft EOS (with the MDI); 
H(M) represents the hard EOS (with the MDI). 


EOS K/ MeV a / MeV B/ MeV o / MeV 
S 200 —356 303 7/6 

SM 200 —390.1 320.3 1.14 

H 380 —124 70.5 2 

HM 380 -129.2 59.4 2.09 


The Coulomb potential can be written as 
follows: 


2 


2 r- 
e 1 zl j 


yuo = i 
Ae (41L) 4L 


1 GQ) 
where 7, = lr, -r| represents the relative distance of 
And L is the so-called Gaussian 
wave-packet width for nucleons which a constant 
L=2.16 fm’ is used here. 

The Yukawa potential can be written as follows: 


two nucleons. 


ay, IDF + exp(Lm?) 


izj fj 
xlexp(mr, Jerfe(dLm-r, / 44L) (4) 
—exp(nr, Jerfe(VLm + 1, /V4L)] 


where V,=0.0024 GeV, m=0.83. 
The symmetry potential can be written as 


follows: 
C l 
oe = Uy -u o; + 2u’6,] 
n (5) 
u 2) 5 2h) 
Py (P); 


where Csym is the strength of symmetry energy at 
saturation density, y is the symmetry energy stiffness 
coefficient, and the symbol “+” for neutrons, “—” for 
protons. 

The MDI can be written as follows: 


al +e ) (6) 


Po 


U = dita Fie 


where 5=1.57 MeV, e=500 Bew . 


2.2 Formulism and Method 


The dipole moment of GDR in coordinator space and 


momentum space can be defined as follows°7”!: 


DRepg O = IR (0) - Ry (7) 


NZ BO Py ©, 
Ah Z N 

R,(@) and R(t) 
Mass (CM) of protons and neutrons in coordinator 


DK pp (‘)= = 


(8) 


Here, represent the Center of 
space, respectively; P,(t) and P(t) represent the 
CM of protons and neutrons in momentum space, 
respectively. 

By the Fourier transformation of the second 
derivative of DR», (tf) in respect of time: 


DR"(@) = ha DR"(t)e'dt (9) 


one can get the photon emission probability for energy 
E, =ho as follows: 


dP 


DR 10 
dE, ~ Bari awe ‘(a (10) 


Similarly, the dipole moment of PDR in 
coordinator space and momentum space can be 
defined as follows: 


NZ 
DRopp (£) =— 


(11) 


DK pop (t) = aa lee s0) (12) 


v 


Here, R,(t) and R, (t) represent the CM of the 
isospin-symmetric core and the valence neutrons in 


coordinator space, respectively; F(t) and P, (¢) 
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represent the CM of the isospin-symmetric core and 
the valence neutrons in momentum space, respectively. 
Then the photon emission probability of PDR can be 
also obtained. 

In PDR calculation, a key point is how to 
choose the valence neutrons in the simulation. In our 
calculations, we choose the neutrons with the farthest 
distance from the CM of all nucleons as the valence 
neutrons in the initial state, and extract their oscillation 
character in the following dynamical process. 

The fraction of Energy-Weighted Sum Rule 
(EWSR) contained in the PDR relative to that located 
in the GDR region can be written as follows: 


Binding Energy / MeV-A* 


RMS / fm 


PDR 
EWSR% = ———"__ 
GDR, + PDR 


m 


x100% (13) 


where the PDR,, is EWSR of PDR and GDR,, is 
GDR, i.e.: 


GDR 5 AEXE (14) 
m = E 
1 Y GDR 
& dP 
PDR,, = >.(——) AExE (15) 
EI ” PDR 


— Event 1 
—— Event 2 
—— Event 3 
—— Event 4 
—— Event 5 


Time / fmc” 


Fig.1 Time evolutions of the root mean square radii and of the binding energy for some samples of stable nuclei of Ni. Soft EOS 


w/ MDI. 


3 Results and discussions 


In the present low energy transport model simulations, 
the stability check for the nucleus is always important. 
To this end, we have to pick up the stable nucleus in 
the initialization. For an example, the time evolution 
of the root mean square radii and the binding energy of 
some stable nuclei are shown in Fig.1. Clearly, both 
the root mean square radii and the binding energy are 
stable with the time evolution before 400 fm-c '. 

Using the above stable nuclei to do Coulomb 
excitation simulation with the IQMD, one can get the 
spectrum of photon emission. Here 197 Au is chosen as 
the target, and only the GDR or PDR spectrum of the 
collective motion of nucleons inside projectile (Ni 
isotopes) is calculated. By a Gaussian fitting to the 


spectra, the peak energy and strength can be extracted. 
Then their sensitivities to impact parameter, incident 
energy and the strength of symmetry energy can be 
also investigated. The sensitivities of these parameters 
to impact parameter have been given in Fig.2. It shows 
that the peak energies and strengths of both GDR and 
PDR do not change much with the impact parameter 
increasing. Considering that the impact parameter 
changes from 16 fm to 26 fm, which makes the 
Coulomb excitation just change a little, this result is 
reasonable. 

Figure 3 shows the incident energy dependence 
of GDR and PDR parameters for “Ni. 

When the incident energy increases, both the 
peak energies and strengths of GDR and PDR decrease. 
With the increasing of incident energy, the projectile 


030502-3 


chinaXiv:202306.00183v1 


ChinaXivA ESAT! 


TAO Cheng et al. / Nuclear Science and Techniques 24 (2013) 030502 


goes through the Coulomb field of target more quickly 
so that the Coulomb excitation becomes smaller, and 
then the projectile gets less excited. 
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Fig.2 Impact parameter dependence of GDR and PDR 
parameters for Ni. In calculations, we use £,,=600 
MeV/nucleon, Csym=32 MeV, y=1, and the Soft EOS w/o MDI. 
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Fig.3 Incident energy dependence of GDR and PDR 
parameters for “Ni. In calculations, we use b=24 fm, y=1, 
Cyym=32 MeV, and the Soft EOS w/o MDI. 

As well known, symmetry energy is very 
important in nuclear physics and astrophysics. In some 
previous studies, PDR has been shown as a useful 
probe to investigate symmetry energy. Here we also 
investigate the relationship between PDR and 
symmetry energy. In our IQMD model, the symmetry 
energy can be written as: 


2/3 C y 
Fg =125( 2] Sel 2) (16) 
Po 2 A 


where the first term is the kinetic energy term and the 


second term is the potential term, y is the symmetry 
energy stiffness coefficient. We study the relationship 
between the two parameters, i.e., Csym and y in the 
symmetry energy formula, and the parameters of PDR. 


Fig.4 shows the strength of symmetry energy Cyym 
dependence of GDR and PDR parameters for “Ni. 
With the Csym increasing, the peak energy of GDR 
increases, while that of PDR decreases and both 
strengths of GDR and PDR decrease. The y 
dependence of GDR and PDR parameters for “Ni is 
shown in Fig.5. With y increasing, the peak energy of 
GDR decreases, while that of PDR increases and both 
strengths of GDR and PDR increase. 
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Fig.4 Csym dependence of GDR and PDR parameters for SN 
In calculations, we use E;,=600 MeV/nucleon, b=24 fm, y=1, 
and the Soft EOS w/o MDI. 
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Fig.5 y dependence of GDR and PDR parameters for Ni. In 
calculations, we use £;,=600 MeV/nucleon, b=24 fm, Cyym=35.2 
MeV, and Soft EOS w/o MDI. 

Figures 4 and 5 show different behaviors of the 
peak energy of GDR and PDR: i.e. the peak energy of 
GDR increases with the Csym increasing, and decreases 
with the y increasing, For PDR, the peak energy 
decreases with the Csym increasing, and increases with 
the y increasing. To understand these questions, Fig.6 
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shows the symmetry energy versus Csym and y. Since 
we use Coulomb excitation method in our calculation, 
the density of nucleon in projectile is mainly in 
sub-saturation area. By the line for guiding eyes in 
Fig.6, one can see that the symmetry energy increases 
with Csym increasing when y is always fixed to be 1 
(Fig.6a), and decreases with the y increasing in 
sub-saturation area (Fig.6b). So the peak energy of 
GDR has a positive correlation with symmetry energy, 
but for PDR an anti-correlation with symmetry energy 
is shown. 
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Fig.6 Esym vs. p/po with various Csym or y. (a) y is fixed to be 1; 
and (b) Csym is fixed to be 35.2 MeV. 
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Fig.7 EWSR vs. the derivative of the symmetry energy at 
saturation (L). 

In previous studies, a good correlation between 
L (the derivative of the symmetry energy at saturation) 
and the fraction of energy-weighted sum rule (in 
percentage) exhausted by the PDR in Ni has been 


[28] [28] 


built-up“. Different mean-field calculation results 


for ®Ni are displayed in the Fig.7, where the solid circles 


represent those results. From the symmetry energy 


Eq.(16), one can get the Z as follows: 
3 
baint VaT (17) 


12 


@ Sn Isotopes 
yz Calculation of Ni Isotopes 
10 |--- Fitting of Calculation 


Fig.8 EWSR(PDR,, /GDR,,% ) vs. N/Z. In calculations, we 


use Ej,=100 MeV/nucleon, b=24 fm, Cyym=32 MeV, y=1, and 
the Soft EOS w/o MDI. 


From the above equation, L increases linearly 
with either Csym or y when another is fixed. Our 
calculation results are shown in open stars, solid stars and 
open diamonds, respectively, for the L deduced from the 
different Csym for y=1 (soft EOS w/o MDI), from the 
different y where C.ym=35.2 MeV (soft EOS w/o MDI), 
and from the different y where Cyym=35.2 MeV (soft 
EOS w/ MDI). It is found that EOS with MDI gets 
higher fraction of EWSR compared with that without 
MDI. This may be due to the effective mass reduced 
by MDI. The curves are their polynomial fits, from 
which one can estimate the range of L and Csym for our 
IQMD model. The box in the figure shows the 
constraints from experiments and other calculations”. 
Within this box Csym and y in IQMD model can be 
limited from the crossing-points of of fitting curves. In 
the case of y=1 for the Soft EOS w/o MDI, one can get 
Cyym€[16.9, 42.9] MeV if we take Le[50.3, 89.4] as 
Ref.[8]. In the case of Csym=35.2 MeV for the Soft 
EOS w/o MDI, we obtain ye[0.61-1.22] and 
Le[57-89.4] MeV. In the case Csym=35.2 MeV for the 
Soft EOS w/ MDI, ye[0.78-1.22] and Le[66-89.4] 
MeV. Based on the above results, the symmetry energy 
with too soft or too stiff density dependence should be 
excluded. 
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Finally, we plot the fraction of EWSR 
( PDR,, /GDR,,% ) versus N/Z for our results 
together with Piekarewicz’s"”! in Fig.8. With the ratio 
of neutrons and protons increasing, our calculation 
results of Ni isotopes show a linear increasing; but for 
Piekarewicz’s results of Sn isotopes within the 
framework of the relativistic random phase 
approximation, it shows a linear increasing firstly, and 


then a slightly linear decreasing. 


4 Conclusion 


In this paper, we studied PDR and GDR of Ni isotopes 
by the Coulomb excitation in the framework of IQMD 
model. By the Gaussian fitting to the photon spectra, 
we obtain parameters of GDR and PDR, such as the 
peak energy and strength. The sensitivities of these 
parameters to impact parameter, incident energy, and 
the symmetry energy are studied. The results show that 
the peak energies of GDR and PDR are correlated with 
symmetry energy and while EWSR% is anti-correlated 
with the derivative of the symmetry energy at 
saturation. Combining with the data and other 
calculation, we constrain the y coefficient and L 
parameter for the IQMD with the different symmetry 
energy parameters, i.e. for the soft EOS w/o MDI, the 
y is in 0.61-1.22 region and L is in 57-89.4 MeV, and 
while for soft EOS with MDI, y is in 0.78-1.22 region 
and L is in 66—89.4 MeV. In the other word, too stiff or 
too soft symmetry energy form is excluded from the 
present calculations. At the last, the N/Z dependence of 
PDR and GDR parameters of Ni isotopes was studied. 
It shows that the EWSR 
(PDR, /GDR,,, % ) increases linearly with the N/Z. 


fraction of 
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